Entanglement of Indistinguishable Particles 



Toshihiko Sasaki 1 . Tsubasa Ichikawa 2 , and Izumi Tsutsui 1 ' 3 
1 Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
2 Research Center for Quantum Computing, Interdisciplinary Graduate School of Science and Engineering, 
Kinki University, 3-4-1 Kowakae, Higashi- Osaka, Osaka 577-8502, Japan 
3 Theory Center, Institute of Particle and Nuclear Studies, 
High Energy Accelerator Research Organization (KEK), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 



o 

O 

Q 

(N 

^— > ■ 

c ■ 
ctf ■ 
=5 



> 



o 
o 



X 



We present a general criterion for entanglement of TV indistinguishable particles decomposed into 
arbitrary s subsystems based on the unambiguous measurability of correlation. Our argument pro- 
vides a unified viewpoint on the entanglement of indistinguishable particles, which is still unsettled 
despite various proposals made mainly for the s = 2 case. Even though entanglement is denned 
only with reference to the measurement setup, we find that the so-called i.i.d. states form a special 
class of bosonic states which are universally separable. 



PACS numbers: 03.65.Ud, 03.67.-a, 03.67.Lx. 



I. INTRODUCTION 



Since its first recognition in the seminal EPR and 
Schrodinger's papers [HQ, quantum entanglement has 
long been seen as the most distinctive trait of quan- 
tum theory. Notably, it underlies nonlocal correlation 
in composite physical systems, invoking various concep- 
tual questions on the foundation of physics and, at the 
same time, offers a key resource for quantum informa- 
tion sciences. In view of this, we find it rather puzzling 
that the very notion of entanglement still eludes a for- 
mal, let alone intuitive, understanding, especially when 
the system admits no apparent decomposition into sub- 
systems. This occurs typically in systems of indistin- 
guishable particles (i.e., fcrmions or bosons) with which 
actual realizations of entanglement - via photons, elec- 
trons or composite particles such as hydrogen atoms - 
have been implemented mostly today. 

To see the nontrivial nature of entanglement, take, for 
example, the familiar TV = 2 particle Bell states, 



|*) = i=(|0) 1 |l} 2 ±|l} 1 |0} 2 ) ) 



(1) 



with |0) fc and |l) fc being the orthonormal qubit states 
of particle k = 1,2. These are prototypical entangled 
states for distinguishable particles, but if the particles 
are indistinguishable, the labels k are no longer usable 
for classifying the measurement outcomes to define the 
correlation. A remedy for this is to consider remotely 
separated particles by introducing the position degrees 
of freedoms directly for ((T|), but this does not yield any 
nontrivial correlation, a property known as the cluster 
separability 0, 13) ■ Clearly, a physically motivated and 
mathematically solid definition of entanglement is needed 
for general composite systems including those of indistin- 
guishable particles. 

Recently, Ghirardi et al. gave a possible definition 
of separability (non-entanglement) for TV particle systems 
based on the criterion that, if one can deduce a complete 
set of physical properties (CSP) pertaining to a subsys- 



tem, then the state is separable with respect to the sub- 
system and the rest. This criterion derives from the de- 
mand that, in a separable state, all physical quantities in 
the subsystem have elements of reality in the EPR sense 
[H . Independently, Zanardi et al. Q presented a criterion 
for uncovering a tensor product structure (TPS) in the 
Hilbert space upon which entanglement can be defined. 
The criterion demands the existence of subalgebras rep- 
resenting the obscrvables of the subsystems, which are 
measurable, independent, and complete to form the en- 
tire set of measurable obscrvables in the system. Yet an- 
other criterion has been proposed by Schliemann et al. [|[ 
and others Q particularly for indistinguishable particles 
using the (Schmidt or Slater) rank of state decomposi- 
tion, which is related to the standard measures of entan- 
glement such as the von Neumann entropy. 

These criteria for (non)cntanglcment arc rather differ- 
ent from each other and, not surprisingly, do not com- 
pletely agree on deciding which states are separable, with 
an example being the TV = 2 bosonic 'independently 
and identically distributed' (i.i.d;) state \cf>} 1 \cf>} 2 (for an 
attempt of reconciliation, see [10|). More recently, the 
present authors furnished a criterion for the decomposi- 
tion of an TV fcrmionic system into s arbitrary subsystems 
[Tlj . where we find that the orthogonal structure intro- 
duced to distinguish the subsystems in Q corresponds 
precisely to the choice of observables with which corre- 
lation is defined. In other words, entanglement can be 
defined only relative to the measurement setup and it is 
highly non- unique Q. Under these circumstances, one 
is naturally led to ask if there is any coherent picture of 
entanglement prevalent among these criteria. 

The purpose of this article is to provide a positive an- 
swer to this. Namely, we show that all these criteria can 
be put into a larger perspective consisting of two descrip- 
tions of the system, one for the measurement outcomes 
and the other for the provisional states of the system. 
The gap between the two descriptions, which lies at the 
root of the apparent disagreement, can be filled by an 
isomorphism between the two, providing a unified view- 
point of entanglement for indistinguishable particles. Un- 
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like the previous analyses, entanglement can be treated 
equally for the fcrmionic and bosonic cases here. We also 
find that the i.i.d. states for general N form a special 
class of bosonic states which are universally separable 
irrespective of the choice of measurement setup. 

II. ENTANGLEMENT IN 
MEASUREMENT-BASED DESCRIPTION 

To define entanglement as an attribute to gener- 
ate nontrivial correlation among subsystems, we first 
need an appropriate set of physical observables associ- 
ated with the subsystems for which the correlation in 
their measurement outcomes can be evaluated unambigu- 
ously. To discuss the situation explicitly, we consider the 
case where the total system breaks into s subsystems 
Ti, . . . ,r s and assume that to each I\ we have a com- 
plete set of commuting observables Ci which are all im- 
plementable in the measurement to determine the state 
of the subsystem. Let Ci be the set of observables (self- 
adjoint operators) containing the set Ci. The collection 
of states of the subsystem I\ describing the measure- 
ment outcomes form a Hilbert space ■H mes (r i ) in which 
Ci is represented irrcducibly. Assuming further that the 
measurements of the observables Ci can be performed in- 
dependently for all i = 1, . . . , s, we find that the set C 
of observables in the total system is given by C = ®iCi. 
Accordingly, the state space of the system describing the 
measurement outcomes is given by the tensor product, 

s 

n mcs = (g)H mcs (r l ). (2) 

i=l 

The TPS of the total space -H mcs in © allows us to define 
the entanglement by the conventional way, that is, if the 
measured state G % mcs admits the product form 

s 

l*>=0W<i>iv IV^r, e H mcs (r,), (3) 

i=l 

then it is separable; if not, it is entangled. Evidently, 
since the separable state ([3]) yields definite outcomes for 
the measurement of observables in a properly chosen set 
Ci in \ipi) r . for all i, it possesses a CSP [5[. 

Meanwhile, in the total space % mes the observable Oi G 
Ci is expressed by 

i— 1 s 

Oi = (g) lj <g> Oi ® (g) lj (4) 

3=1 j=i+l 

where lj is the identity operator in H mes (Tj). The afore- 
mentioned independence is then assured trivially by 

[d i} dj]=0 for i^j. (5) 

The observable O G C corresponding to the simultaneous 
measurement for the subsystems is then given by 

i=i i=i 



Denoting the set of such operators by 7" mcs c C, we see 
that any O G 7~ mcs has a factorized expectation value for 
the separable state |^) in ([3]): 

s 

(*|d|*) = n<*|6il*>- (7) 

i=l 

The properties ([5]) and ©, together with the implc- 
mentability assumption, constitute the formal conditions 
to realize a TPS in Q- Note that, in our measurement- 
based description, the TPS appears as a direct conse- 
quence of the construction. 

III. ENTANGLEMENT IN PROVISION-BASED 
DESCRIPTION 

The entanglement in the measurement-based descrip- 
tion is related with the measurement outcomes directly, 
but the conventional treatment of indistinguishable par- 
ticles employ the framework of the provisional Hilbert 
space of distinguishable particles for the description of 
states with appropriate restriction required by the statis- 
tics of the particles. Here, the description is not directly 
related to the measurement outcomes, and the restricted 
space of states does not admit a TPS in any obvious man- 
ner. In physical terms, the measurement outcomes of ob- 
servables, such as spin, cannot be attributed to those of 
a particular particle due to the indistinguishability, and 
the formal structure of the state fails to signify the cor- 
relation as exemplified by flU). To fill the gap, we need a 
prescription to bridge the two descriptions. 

For definitcness, let us label the N particles by the 
integer set N = {1.2, ••• . N}. Each of the particles is 
characterized by an n- level state, i.e., the state space of 
the fc-th particle is Hk — C n for all k. Let {|ej)} be a 
complete orthonormal basis in C". By the isomorphism 
among the constituent spaces Hk, any pure state |\&) in 
the provisional space T-L = ®fcejV ^ fc °f * ne total system 
can be written as 

N 

l*)= E ® K>fc > ( g ) 

i\ ,22 ' An k—i 

where ^t 1 t 2 ...i JV G C and {|ei fc ) fe } is the complete or- 
thonormal basis in Hk isomorphic to {|e,}}. 

To incorporate the indistinguishability of the particles, 
consider an element a G &n of the symmetric group &n 
associated with the permutation k — > a(k). In the 
element is represented by a self-adjoint operator ir a with 

N 

i\ ,22 r" An k—1 

From 7T CT , both the symmctrizcr and the antisymmctrizer 
are defined as 

s = ^Y. iT °> A= ^J2 s ^w, (io) 
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where sgn(cr) is the signature of the permutation a . The 
Hilbcrt space of the total system of N bosons (fermions) 
is the subspace of H consisting of symmetric (antisym- 
metric) states. Putting X = S for bosons and X = A for 
fermions, they are 



H x = [H] x ~{X\V)\\*)eH} 



(11) 



To introduce the decomposition into subsystems in the 
total system, we consider a partition T of the integer set 
Af into non-empty and mutually exclusive sets I\ C J\f, 

s 

r = {r*K =1 , Ur i= ^, r 4 nr 3 = for;^j.(i2) 

i=l 

For specifying the subsystems of indistinguishable parti- 
cles, only the cardinality |I\| of I\ matters. Note that 
there is no apparent TPS in Ji x with respect to T, and we 
need to somehow find an embedding of the measurement- 
based description in the provision-based description. 

This embedding is handled usually by considering the 
positions of particles to gain a fictitious distinction be- 
tween the particles. For the distinction to be unambigu- 
ous, the measurements of the subsystems should be per- 
formed remotely from each other, and this amounts to in- 
troducing an orthogonal decomposition in the 1-particle 
Hilbcrt space (after accommodating the position degrees 
of freedoms). More generally, the embedding requires 
an orthogonal structure V which is a set of subspaces 
Vi C C™ mutually orthogonal to each other with respect 
to the innerproduct of C™, 

V = {V i \U, V i ±V j faH^j. (13) 

Together with the orthogonal complement, 

V a = (Vi 8 V 2 © • • • © Vi) X , (14) 

the set V furnishes an orthogonal decomposition of C™. 
The physical idea behind this is that these orthogonal 
spaces correspond to mutually independent measurement 
of subsystems such that, given a measurement setup, the 
subsystem Ti is susceptible only for the measurement of 
particles k £ Tj residing in Vi . If we denote the subspace 
Vi in Hk by Vi(Hk) C Hk, then the actual Hilbert space 
describing the measurement outcomes for Ti is given by 



n x {T t ,Vi) 



Vi(H k ) 



(15) 



x 



Here, as in (fTTj) . the index X attached to [*] implies 
that it is the subspace of * invariant under X associ- 
ated with the symmetry group &k with K being the 
cardinality of *, and in (TT5")) we have K = |r.;|. Clearly, 
Hx(Fi,Vi) is the actual space of states determined from 
the measurement and, therefore, corresponds to - H mes (r i ) 
in the measurement-based description @ where the state 
\ipi)r- G HxiThVi) is identified with |^) r . € U mes (Ti). 



From the description for the subsystems, we obtain the 
Hilbert space of the total system by 



n x (r,v) = 



(16) 



x 



Note that, due to the (anti)symmetrization X, the re- 
sultant space T-L X (T,V) has no TPS with respect to 
the decomposition T and, hence, no obvious correspon- 
dence with % mcs in In spite of this, the two 
spaces can be made isomorphic based on the identifica- 
tion Hxi^i, Vi) = T-L mes (Ti) mentioned above. Indeed, it 
is attained, with this identification, by the map, 



fx ■ n r 



defined by 



\Ux{Ti,Vi)^U x {T,V), (17) 




MX\ 



i/Ti 



(18) 



i=i 



with the normalization factor M := N\/Y[i\^i\- Ob- 
viously, the map f x is surjective by construction, and 
to see the injectivity, we note that, thanks to the or- 
thogonal structure V in (|13[) . the innerproduct is invari- 

i \^i)r t II = 



12(. It follows that 



ant under the map 

\\\/~MX ®i |V't)r- Hi which ensures the injectivity of the 
map and hence the isomorphism (see FigfT]). 

The isomorphism (fTTj) induces the correspondence be- 
tween the observables in the two spaces. If Oi are the ob- 
servables in H^(r,, Vi) for i = 1, . . . , s, then the observ- 
able in H X (T, V) for simultaneous measurement reads 



O = MX ( (g) Oi J X. 



(19) 



The set of all such operators defines a subset T(r, V) 
of observables in T-L X {T,V). With the identification of 
the observables Oi between H x (Ti, Vi) and "H mes (rj), the 
induced isomorphism for the observables corresponding 
to (fTHJ) is given (by abusing the symbol) by 




Oi = MX 




(20) 



This also implies the isomorphism between T™ 03 and 
T(r, V) through the correspondence 0^0. 

In the provision-based description, the criterion on the 
entanglement of indistinguishable particles then emerges 
as follows. Given an arbitrary (normalized) state |^) € 
T-L x , we first decompose it as 



|*) = |*(r,y)) + |vi/(r,y) ± ), 

according to the orthogonal decomposition, 
Hx=H x {r,V)® Hx(T,V) x . 



(21) 



(22) 
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FIG. 1: A diagrammatical representation of the spaces men- 
tioned in (|17l) and 1|22|). In the total space Hx, we find the 
subspace Hx(T, V) isomorphic to ®iWx(yi, Vi) which has a 
TPS. The latter is then identified with the space H mcs de- 
scribing the measurement outcomes. 



Since (^(r, V)^) has a vanishing support for the observ- 
ables in 7~(r, V) and is filtered out by the measurement, 
the only part significant for correlation is |W(r, V)). 
Thus, for studying correlations in the measurement out- 
comes ignoring the events which are not detected in the 
setup, one uses the remormalizcd state || ^(r, V)) \\ = 1. 
We now see that, if the observable part takes the form, 



\*<T,V)) 



(23) 



then for any O € T(r, V) we have the factorization, 



(vi/|oi*)=n<*io l i*> 



(24) 



in analogy with ([7]). Since the converse is also true, we 
learn that the state |^) is separable if and only if the 
part |^(r, V)) in (j2Tj) admits the (anti) symmetrized di- 
rect product form (|23| : if not, it is entangled. In more 
simple terms, to examine the separability of a given state 
| 'J) € 7ix, we just concentrate on the observable part 
^(r, V)) and then strip it off the projection X (and per- 
form necessary renormalization) to obtain, via the identi- 
fication in Fig (TJ the corresponding state |^E' m08 (r, V)) <G 
H mcs describing the measurement outcomes directly. In 
the case (J23J) we find |* mos (r, T^)) = ®| =1 \tpi) T ., which 
is factorizable and hence separable. As is evident from 
the explicit dependence on V, the entanglement of the 
state is determined only relatively with respect to the 
measurement setup. 

Generalization of our argument to mixed states is 
straightforward. Given a density matrix p on Hx, one 



o 



0.5 



2it 



(measurement setup parameter) e 



FIG. 2: The concurrence C of |* mcs (r, V)} as a function of 
the angle 8. The variation shows that the state undergoes the 
change from complete separation to maxmial entanglement 
by the adjustment of the setup angle 9. 



can decompose it as 



p(T,V) 



(25) 



where p(T, V) and p(T, V) 1 - are (unnormalizcd) density 
matrices on Hx(F, V) and Hx(F, V) 1 - , respectively. By 
virture of the isomorphism fx, the separability crite- 
rion for the mixed distinguishable systems [l3j can be 
utilized for the density matrix p mes (T, V) which is de- 
fined from f^. 1 (p(I\ V)) with a suitable rescaling to fulfill 
Tr p mes (r, V) = 1. We then find that a mixed state p on 
"Hat is separable under our measurement setup specified 
by r and V if p mos (r, V) admits the form, 



3 (r,v) 



Pa 



I*: 



(26) 



where l^™ 08 ) € f-[ mes are separable pure states and {p a } 
satisfies p a — 1 and p a > 0. 

It should be clear by now that since a state of indis- 
tinguishable particles, either it is pure or mixed, can be 
mapped to a state in H mcs , the entanglement of the state 
can be studied in terms of the standard entanglement 
measures developed for distinguishable particles. This 
will be demonstrated next. 



IV. MEASUREMENT SETUP DEPENDENCE 
OF ENTANGLEMENT: EXAMPLES 

In order to evaluate explicitly the dependence of en- 
tanglement on the measurement setup we choose, let us 
consider the case of iV" = 2 fermions, each possessing 
n = 4 dimensional constituent space, given in the state, 



|e 3 ) 2 e?U 



(27) 



We wish to examine if, and to what extent, the state 
is entangled under the partition V = {{1},{2}} when 
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our setup is 'rotated' among the set of basis {|ei) , |e 4 )} 
and {|e2) , |e3)}. To this end, we adopt the orthogonal 
decomposition (which defines the measurement setup) 
V = {Fi.Va} with 

Vx = span{|ei) , \e' 2 )}, V 2 = span{|e 3 ) , \e' A )}, (28) 

where 



\) =c |ei) + s |e 4 ) , \e' 2 ) = -s |e 3 ) + c |e 2 ) 
i) = c |e 3 ) + s \e 2 ) , |e 4 ) = -s |ei) + c |e 4 ) 



(29) 



and we have used the shorthand c = cos 6*, ,s = sm9 
to express the rotation with angle 9. According to the 
decomposition (|21[) . the measurable part turns out to be 

(30) 



|*(r,y)) = ^(c 2 |e' 1 ) 1 | e3 ) 2 - s 2 | e ^ 1 K> : 



up to a constant. We then map it to the corresponding 
normalized state in % mes as 

1 



|* meB (r,v~)) = 



Vc 4 + s 4 



(c 2 le 



l)l \ e 3)l 



s 2 \e' 



2/1 1^4/2) 



(31) 

The amount of entanglement may be evaluated by the 
(squared) concurrence, 



c(|* mes (r,i/))) = 2 



1 - Th (Tr 2 |* r 
4 



(\f mes |)" 



(tan 2 9 + cot 2 



(32) 



and the result is depicted in Fig. [2l We find that the 
state (|27|) is strictly separable at 9 — nix /2 and maximally 
entangled at 9 = (n + l/2)7r/2 for integer n, and it can 
take any intermediate values C in between. 

As a second example, we consider the case of N = 2 
bosons with n = 6 prepared in the state, 



I*) 



1 



6 



V6^ 

V 2—1 



(33) 



As before, we study the entanglement for the partition 
r = {{1}, {2}} when the 'rotated' family of entanglement 
setups are considered, which are now provided by 



Vi = span{|ei) , \e 2 ) , |e 3 )}, 
with 

\<& = U\ei 



V 2 = span{|e 4 ) , \e 5 ) , leg)}, 
(34) 



U e U(6). 



(35) 



Analogously to Eq. (j31j) . we can find the corresponding- 
state |f mcs (r, V0) in n mcs . This time, however, instead 
of simply evaluating the concurrence we study the extent 
of variation in the state |^E' mos (r,y)) that can arise un- 
der distinct measurement setups obtained by altering the 
unitary matrix U in (|35|) . To do this, wc first implement 
the Schmidt decomposition for the state |\l/ mes (r, V)) as 



|* mos (r,^)) 



i=l 



i+3/3 



>3> 



3 

E 



A 2 = l, (36) 




0.333 



FIG. 3: The distribution of the Schmidt coefficients of the 
state (|36p with the ordering Ai < A2 < A3 (A3 is not shown 
because it can be determined from other two) obtained under 
random setups of measurement for the state (|33[) provided by 
10 6 distinct unitaries U. The diagonal lines represent Ai = A2 
and A2 = A3, whose intersection corresponds to the maximally 
entangled state. The points fill out basically the entire region 
of the triangle, although in our random generation the density 
becomes scarcer for states which are highly entangled. 



where {|e")}f =1 , {|e-')}f =4 are the Schmidt bases 
{l e i')}f=i: {l e i')li=4 eac h defined within the measurable 
subspaces Vi, V 2 . We then observe the distribution of 
states which is invariant under local unitary operations 
from the distribution of Schmidt coefficients. FigOshows 
the values of for i = 1, 2 obtained by a random gener- 
ation of U , which suggests that by tuning U properly the 
state (|33j) can furnish virtually any possible states which 
are discriminable by the Schmidt coefficients. 

Despite the relative nature of entanglement with re- 
spect to the measurement setup, there exists a special 
class of states in the bosonic case which are separable 
under all measurement setups. These are the i.i.d. pure 
states IvP) € Hs defined by 



A? 



I*)=®l 



Ik ' 



(37) 



fc=l 



To see the universal separability of the state, we decom- 
pose \cj)) k according to Eqs. fTS")) and (Tl4]) as 

s 

l<^ = El^>fc< \<Pi) k eVi(H k ). (38) 

i=0 

Plugging this into Eq.([37|, we obtain Eq.([2l]) with 



(39) 



\*(T,V)) = VMSQ$\1> i ) rt , 
i=i 



where 



\i>i 



r, 



oc 



(40) 



feeri 
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Since the part V)), if no n- vanishing, belongs to the 

class Eq. (|23|) . we see that the i.i.d. states (|37|) are separa- 
ble. Further, since this is true for any choice of (T, V"), the 
separability holds irrespective of the measurement setup. 
Interestingly, one finds that, for N = 2, n > 4, the con- 
verse is also true: states which are universally separable 
must be the i.i.d. states. 



V. CONCLUSION AND DISCUSSIONS 

In summary, we have presented a general criterion for 
entanglement of an indistinguishable N particle system 
decomposed into s subsystems based on the unambiguous 
measurability of correlation. The point is that, although 
the Hilbert space Hx of the system does not admit a 
TPS, one can find a subspace %x(P, V) C T-Lx which has 
a TPS and is directly related to the space % mes describing 
the measurement outcomes. Since % mes has a common 
structure with the space of distinguishable particles, our 
approach allows us to treat indistinguishable particles on 
the equal basis with distinguishable ones. Consequently, 
the handling of states without considering the effect of 
(anti)symmctrization practiced regularly in quantum op- 
tics is found to be safe as long as it deals with the space 

^mes 

The structure of fi mes also implies that the standard 
measures of entanglement devised so far can be used 
equally for the indistinguishable case; for instance, the 
monotonicity of entanglement measures with respect to 
local operations and classical communications (LOCC) is 
preserved under the mapping fx- This is shown by ob- 
serving that all ingredients of LOCC for distinguishable 
particles Q have their counterparts in Wx(r, V) pro- 
vided by the application of fx- Since generalized mea- 
surements, POVM, can be built from some of the ingredi- 
ents of LOCC (Naimark's theorem Q), the mapping fx 
induces the analogues of generalized local measurements 
\nUx{T,V). 



As stated in the introduction, for bosonic systems 
the characterization of separability has not been done 
uniquely in the literature, and in fact it is mentioned in 
(Theorem 3.6) that there are states which are sepa- 
rable from the criterion of 0, Q but entangled from that 
of [||. We observe, however, that the theorem used is 
based on a unitary transformation between two different 
bases of the one-particle state space, which amounts to a 
change in the orthogonal structure V in our language. As 
explicitly shown in the previous section, such a change 
gives rise to states with different amounts of entangle- 
ment observed in the altered settings, and in this sense 
one can regard the apparent discrepancy as just a reflec- 
tion of the relative nature of entanglement. 

Finally, we mention that our approach can also be 
applied, rather trivially, to a system consisting of both 
bosons and fermions. Indeed, we may first treat bosons 
and fermions separately in the provisional spaces Jig and 
Ha in our approach, and then combine them together to 
form the total provisional space Hs <£> Ha- Since Hs 
and "Ha have their own isomorphism fs and /a defined 
from the orthogonal structures equipped with them, it is 
evident that through the combination of f^ 1 and /7 1 
we obtain the corresponding state space which has a 
TPS and can be identified with the total measurement 
space H mcs . This reasoning can be extended in principle 
to more complex systems consisting of several distinct 
species of fermions and bosons, e.g., those describing in- 
teractions between matter and gauge mediators such as 
photons and gluons. 
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